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ABSTRACT 

Recent observations of the ground state transition of HDO at 464 GHz towards 
O . the protoplanetary disk of DM Tau have detected the presence of water vapor 

^ . in the regions just above the outer disk midplane (Ceccarelh et ah 2005). In 

. the absence of non-thermal desorption processes, water should be almost entirely 

Q . frozen onto the grain mantles and HDO undetectable. In this Letter we present 

^ , a chemical model that explores the possibility that the icy mantles are photo- 

^ . desorbed by FUV (6eV < hu < 13.6 eV) photons. We show that the average 

y—i '. Interstellar FUV field is enough to create a layer of water vapor above the disk 

! midplane over the entire disk. Assuming a photo-desorption yield of 10"^, the 

'. water abundance in this layer is predicted to be ~ 3 x 10~^ and the average H2O 

^ ! column density is ~ 1.6 x 10^^ cm~^. The predictions are very weakly dependent 

^ ! on the details of the model, like the incident FUV radiation field, and the gas 

! density in the disk. Based on this model, we predict a gaseous HDO/H2O ratio 

^ ! in DM Tau of ~ 1%. In addition, we predict the ground state transition of water 

at 557 GHz to be undetectable with ODIN and/or HSO-HIFI. 

o ■ 

^ I Subject headings: Stars: formation - Stars: protoplanetary disks - Stars: Pre- 

■ main-sequence - ISM: molecules - 



Introduction 



Water is a key ingredient in many astrophysical environments, including circumstellar 
disks. It is the main constituent of the icy mantles coating interstellar grains, and can 
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be one of the most abundant molecules in the gas phase. For these reasons, water may 
dominate both the chemistry (being a major oxygen reservoir) and the thermal balance (i.e. 
cooling) of the gas component of protostellar disks. In addition, water vapor and water ice in 
protoplanetary disks are the major reservoirs of water later to be found in planets, asteroids 
and comets and, therefore, are of great relevance for understanding the origin of the Solar 
System and the distribution of volatiles within it. Unfortunately, the observation of water 
molecules in interstellar space is nearly impossible from the ground because of the absorption 
by the terrestrial atmosphere. So far, just a handful of space-based instruments (ISO, SWAS 
and ODIN) have carried out observations of water vapor lines in interstellar molecular clouds, 
but these instruments did not have the sensitivity to observe water protoplanetary disks. 
However, the up-coming ESA/NASA mission Hcrschel will provide a sensitive instrument, 
HIFI, for these studies (Dominik & Ceccarelli 2005). Also, SOFIA has potential for water 
observations. 

Meanwhile, HDO, as the most important isotope of water, is currently the best ground- 
based probe of the presence of water in astrophysical environments. Recently, Ceccarelli 
et al. (2005) reported the first discovery of HDO (and therefore also the first discovery of 
water) in a protoplanetary disk, namely DM Tau. The observations show an absorption line 
of gas phase HDO with a line-to-continuum ratio of 0.9 which translates into an average 
column density through the disk of ~ 1.6 x 10^^ cm~^. In the region where HDO is located, 
its fractional abundance is estimated to be about 3 x 10~^. If the abundance of gaseous H2O 
is similar to the value observed in molecular clouds and protostars, ~ 10"^, HDO/H2O ratio 
is 0.01. The fact that the line is a deep absorption fine immediately indicates that HDO 
must be present above the midplane (from where the sub-millimeter continuum is emitted). 
It must also be present especially in the outer radial regions (> 600 AU) of the disk since the 
continuum originates from an extended region. The presence of a large column of HDO in 
this location was a surprise, since the grain temperatures in this region arc below 25 K, and 
H2O and HDO should be almost entirely frozen out onto grains. Chemical models including 
thermal desorption of water ice and gas- phase formation of H2O (Aikawa et al. 2002; van 
Zadelhoff et al. 2003) predict a H2O column of 10^^ cm~^ at about 400 AU, only a factor of 5 
larger than the observed HDO column. Willacy & Langer (2000) included photo-desorption 
by a strong stellar UV field and found H2O column densities of about 5 x lO^^cm"^ at 700 AU. 
An extended component of HDO (and consequently water) therefore indicates a desorption 
agent from grain surfaces. Important possibilities are X-rays originating from the star, and 
penetrating throughout the entire outer disk, and UV photons, either from the star or simply 
due to the interstellar radiation field. 

High energy photons can in principle act to remove H2O molecules from the ice layers 
on dust grains. X-rays are energetic enough to remove a number of molecules, if sufficient 



-3- 



energy can be concentrated close to the surface of the ice mantle on a grain. This idea has 
been explored by Najita et al. (2001) who considered X-ray heating of small grains, or spot 
heating of larger grains. Small grains, or small thermally insulated spots on large grain can 
be heated strongly enough by a single X-ray absorbed photon to lead to thermal evaporation 
of part of the ice mantle. 

An alternative mechanism is desorption by FUV photons (Willacy & Langer 2000). 
An FUV photon absorbed in the surface layer of an ice mantle can lead to the release of 
a molecule into the gas. In this Letter we explore the effect of just the interstellar FUV 
radiation field irradiating the disk surface of DM Tau. We show that such a scenario leads 
to a layer of water vapor in the disk surface to an FUV optical depth around unity. The 
column densities reached in steady state appear to be a viable explanation for the HDO 
absorption line seen in DM Tau. To show this we developed a simple model of the chemistry 
leading to the presence of water vapor in irradiated gas containing grains covered by layers 
of water ice. A much more detailed discussion in the context of a full Photo-Dissociation 
Region (PDR) model, as well as analytical formulas to estimate H2O column densities under 
such conditions will be given by HoUenbach et al (2005, in preparation). Here we report 
the application of this model to the case of DM Tau. A forthcoming paper will present the 
results of a larger parameter space study, which will explore the dependence of the model 
predictions especially on the characteristics of the grains (Ceccarelli et al. in preparation). 



2. Model description 

The model computes the H2O abundance across a protoplanetary disk, as a function of 
the radius and height. For the physical structure we used the dust density and temperature 
profile which fits the Spectral Energy Distribution (SED) of DM Tau (Ceccarelli et al. 2005), 
and shown in Fig. 1. The gas is assumed to be fully coupled with the dust in terms of density 
distribution and temperature. 

Due to the low temperature, the dust grains are covered with an ice layer. In our 
model, UV photons entering the disk surface are absorbed by dust grains, and cause photo- 
desorption of water molecules from the ice. Only UV photons absorbed directly at the 
surface of the ice mantle can desorb a molecule, so that only a small fraction of absorbed 
photons will cause the ejection of H2O. Laboratory experiments at Lyman a wavelength have 
shown that the involved yields are typically between 10~^ and 10~^ molecules per photon 
(Westley et al. 1995b). The yield increases to some extent depending on the UV dose in 
the experiment, apparently because radical formation on the ice surface aids the desorption 
process (Westley et al. 1995a). It is unclear how important this effect is in astrophysical 
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environments, and we adopt a conservative value for the yield of 10~^ molecules per incident 
FUV (6< hu <13.6 eV) photon. The yield probably introduce a factor of a few uncertainty 
into the computations (see below). 

The desorption process at a given location will also depend on the attenuated FUV flux, 
and the rate per unit volume can be written as: 

where Go is the FUV field in Habing (1968) units, /o = 10^ photons/cm^/s is the FUV 
flux for the standard Habing interstellar field (Go—l). Y is the photo-desorption yield. 
A^uv is H2 column density that gives = 1 between 6 and 13.6 eV; we adopted a value 
equal to 1.8 x 10^^ cm'^ (Tielens & Hollenbach 1985). Cgr is the grain geometrical cross 
section, equal to vr a'^^., where we assumed an average grain size of 0.1 fim. rigr is the grain 
number density. The product Ugrngr is approximatively 10~^^ cm~^ times the gas density 
uhj, obtained assuming a gas to dust ratio in mass equal to 1:100. After being released into 
the gas phase, water either freezes out back onto a grain, or it is photo-dissociated by the 
FUV photons at the rates: 

^freeze — 'S'grCgrHgrnHaO < l^th > (2) 

and 

A;phd = GqIq exp(-A^(i?2)/-/Vuv)nH20 (3) 

where S^y. is the sticking coefficient, for which we use 1 (Burke & Hollenbach 1983). < Vth > 
is H2O molecule thermal velocity, /q = 5.1 x 10~-^'^s"-^ (Le Teuff et al. 2000) is the rate of 
FUV dissociation in unshielded Gq — 1 field. 

Gaseous atomic oxygen can freeze out onto grains, where it is assumed to form water 
ice by reactions with hydrogen atoms on the grain surfaces. O atoms may also form water 
in the gas phase via the standard sequence of reactions started by the reaction between O 
and H^. The gas phase formation of water thus proceeds at a rate: 

^gas = hoxrarionjj+ (4) 

where /cform is equal to 8 x 10~^° s~^ cm~"^ (the rate coefficient for the reaction O + H3 — > 
0H+ + H2) times 0.33 (the last factor accounts for the fraction of H3O+ recombinations 
with electrons forming H2O (Le Teuff et al. 2000)). The atomic oxygen abundance is also 
computed by the steady state equilibrium between formation and destruction processes: 
photo-dissociation of gas phase water (Eq. (3) ) for the formation, and formation of gaseous 
water via reaction with H3 (Eq. (4) ) and freezing onto the grains for the destruction of O. 
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We assume that these reactions are the dominating processes, and that all the oxygen not 
contained in CO or silicates is contained in O, H2O, and H2O ice. 

Finally, the H3 abundance is computed as follows. H3 is formed by cosmic ray ionization 
of H2 to form H2", followed by the reaction with H2 to form H3 . In the regions where 
CO molecules are not frozen onto the grain mantles, H3 is destroyed by the reaction with 
CO. Elsewhere, the H3 abundance is computed following the model described in Ceccarelli 
& Dominik (2005), which takes into account all three deuterated forms of Hj^ and solves 
the chemical composition by considering the reactions between H+, e~, grains, and 
isotopologues. 

3. Discussion 

Figure 2 shows the H2O abundance (with respect to H2) across the disk, for a standard 
interstellar UV field (Gq = !)• Figure 3 shows a vertical cross-section of the chemical species 
involved in water formation at a radius of 700 AU. Despite the relatively large densities in 
the disk, water vapor has an abundance of ~ 3 x 10^'' in a large fraction of the outer disk, in 
the layers just above the midplane, where the Ay to the disk surface is lower than ~ 5 mag. 
By midplane here we mean the region where more than 2/3 of CO molecules are frozen onto 
the grain mantles. This occurs at an height of about 230 AU at a radius of 700 AU, and 100 
AU at a radius of 400 AU. Above this height, CO is desorbed thermally^. The gas phase 
water abundance peaks near the surface where the FUV desorbing fiux is high and where 
the abundance is high. In practice, the freeze-out rate of the H2O molecules is larger 
than the FUV photo-desorption rate of the grain mantles only at densities larger than about 
10^ cm~^. At lower densities, assuming water is in the grain surfaces, water vapor is formed 
mostly through direct photo-desorption from mantles, and, to a lesser extent, by gas phase 
reactions occurring among atomic oxygen (a product of H2O photo-dissociation) and H3". 
In the upper layers, atomic oxygen is predicted to be very abundant. The resulting average 
gas phase H2O column density of a face-on disk is 1.6 x 10^^ cm~^. Figure 4 shows how the 
column density is distributed as function of the radius. The column density is remarkably 
constant over most of the disk, although the column does rise slightly with radius. 

The first conclusion of our study is that vapor water can indeed he found with relatively 
high abundances (~ 10~^ — 10~^j and column densities (~ 1.6 x 10^^ cmT'^) in the protoplan- 



^In our model CO molecules are not photo-desorbed from the mantles, although this is plausible, depend- 
ing on the structure of the mantles and where on the grain the frozen CO is located. However, the gas phase 
CO abundance has only a minor influence on the gas phase H2O abundance in our model. 
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etary disks which surround low luminosity protostars illuminated by the interstellar radiation 
field (ISRF), like in the case of DM Tan. This is caused by the photo-desorption of the icy 
grain mantles by the average ISRF. Increasing the FUV field by a factor of ten and hundred 
leads to H2O column densities larger by only factor 1.5 and 2 respectively. Decreasing the 
FUV field by a factor ten results in decreasing the II2O column density by a factor a bit 
more than a factor 2. Therefore, the H2O column density is not sensitive to the addition of 
a possible FUV field from the central source or a nearby hot star (see HoUenbach et al. 2005 
for details and an analytical proof of the insensitivity to Go). The fundamental reason for 
this insensitivity is that in higher fields, the increased photo-desorption of the ice is balanced 
by the increased photo-dissociation of gas phase H2O. 

As discussed above, the value of the photo-desorption yield, y, is only constrained 
within a factor roughly three to ten by laboratory experiments. This uncertainty causes an 
uncertainty in the predicted H2O column density by the same factor. Another important, 
and poorly known parameter that enters in these computations is the dust-to-gas ratio. In 
the standard case we assumed the canonical mass dust-to-gas ratio equal to 1%. If this 
value is increased by a factor ten (which means decreasing by a factor 10 the gas density in 
the plot of Fig. 1), the II2O column density increases by a factor ~ 1.2. On the contrary, a 
decrease by a factor ten of the dust-to- gas ratio leads to a decrease by a factor 1.2 of the H2O 
column density. This is due to the fact that the UV optical depth is determined by the dust 
distribution only and independent of the gas density. As long as photo-desorption, photo- 
dissociation and freeze-out are the dominant processes, the derived densities and column 
densities (as opposed to the abundances) of H2O are independent of the gas density. The 
small changes are due to regions in which the gas formation rout of H2O is important. 
Therefore, also a drastic change in the dust-to-gas ratio does not cause much of a change 
in the H2O vapor column density. Finally, the results will also depend on the grain size 
distribution. If much smaller grains than 0.1/i dominate the grain surface area, the balance 
between the different formation and destruction rates will be shifted. This question will be 
discussed in a foUowup paper. 

In conclusion, we predict a H2O column density in protostellar disks similar to that 
surrounding DM Tau of a few times 10^^ cm~^, rather insensitive on the external FUV field 
and/or the dust-to- gas ratio in the disk. 

Ceccarelli et al. (2005) reported the detection of HDO in DM Tau, with an observed 
column density of ~ 1.6 x 10^^ cm~^. Based on the present theoretical model, the H2O 
column density is predicted to be ~ 2 x 10^^ cm~^, within about a factor of a few, given the 
uncertainties in the yield, FUV field and the dust-to-gas ratio. This imphes a HDO/II2O ratio 
equal to about 1%. This is consistent with estimates of the HDO/H2O ratio in embedded 
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low mass protostars. Parise et al. (2005) measured a HDO/H2O ratio equal to 3% in the 
sublimated ices surrounding IRAS16293-2422, and a much lower (<0.2%) value in the outer 
envelope where the gas phase chemistry dominates the HDO and H2O formation. Note that 
the large HDO/H2O ratio of 3% refers only to the sublimated ices, which are responsible 
for the water abundance 3 x 10~^ in this object (Ceccarelli et al. 2000). Very likely, the 
sublimated ices are only the "last coating" of the mantles where the molecular deuteration 
is the largest, due to the history of ice formation. In agreement with this interpretation, 
NIR observations of HDO ice have, so far, failed to detect it, at a level of HD0/H20< 2% 
(Parise et al. 2003) Therefore, only a small fraction of elemental deuterium is iced on the 
mantles. The HDO/H2O ratio in DM Tau, equal to ~ 1%, is just shghtly lower than the 
value measured in IRAS16293-2422 and attributed to sublimated ices. If the values were 
indeed different, this might point to reprocessing of the ices in the disk - however, for this 
an actual measurement of the vapor water in DM Tau is required. The measurement of the 
H2O column density would also allow to measure the photo-desorption yield. In addition, 
measures of the HDO/H2O ratio in protoplanetary disks will probe the history of the ice 
formation, with important implications for the chemical and isotopic composition of Solar 
System objects. Based on our model, we predict the ground state transition of H2O at 557 
GHz to be in absorption similar to the HDO ground transition at 464 GHz. The line is 
predicted to be highly optically thick, with r ~ 500. Unfortunately, neither the currently 
available ODIN nor the planned HSO satellites will be able to detect the continuum of DM 
Tau, because of the large beam dilution of the signal from the disk. 



4. Conclusions 

The results reported in this Letter lead to the following important conclusions: 

1. An external UV radiation field as weak as the average Habing field (Go=l) can produce a 
layer of H2O on top of a protoplanetary disk (like that surrounding DM Tau) with a column 
of a few times lO^^cm"^. 

2. This column density is insensitive to external parameters like the strength of the radiation 
field: Variations in Gq by two orders of magnitude only change N(H20) by a factor of 2. 

3. Similarly, changing the dust-to-gas ratio (i.e. the gas density at a given dust density) has 
only weak influence on the derived column densities. 

4. Because of the weak dependence on parameters, the measurement of the H2O column 
density in protoplanetary disks is an excellent tool to measure the photo-desorption yield Y. 

5. Combined with previous observations of HDO (Ceccarelli et al. 2005), we predict an 
HDO/H2O ratio in the outer disk of DM Tau of 0.01, with an uncertainty of about a factor 
of a few. 
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6. We predict that the H2O ground state hne at 557 GHz is undetectable in DM Tau with 
ODIN and/or HSO-HIFI observations. 

Finally, we remark that other molecules trapped in the water ices could also be released into 
the gas phase, keeping "alive" the gas chemistry in the layers above the midplane of the 
outer disk. 
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Fig. 1. — The logarithm of the gas density (upper panel) and the temperature (lower panel) 
of the disk surrounding DM Tau, as derived by modeling of the SED (Ceccarelli et al. 2005), 
and assuming that dust and gas are coupled with Pdust/Pgas — 0.01. 
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Fig. 2. — n(H20)/n(H2) as a function of position in the disk. The sohd hne shows where 
the UV optical depth measured vertically to the disk plane is unity: T^y — 1 (see text). 
The dashed lines show iso density contours: from the top to the bottom they are 10^, 10^, 
10^ and 10^ cm~^. The white region was excluded from the calculation because of the low 
density there. 
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Fig. 3. — The abundances and the visual extinction (from the disk surface) as a function 
of the height above midplane at 700 AU from the star. The sohd and dashed hues show 
the abundance of H2O and O, respectively. The dashed-dotted and the dotted-dashed lines 
show the abundance of CO and H;^. The dotted line indicates the Ay- The most important 
contribution to the H2O column density is produced between Ay ~ 1 and ~ 4 mag, and 
a height between 100 and 250 AU. CO freezes out at a height of ~ 230 AU. This causes 
a small increase of the abundance, since CO (along with O) is the main destroyer of 
H3 . The high abundance of CO at very low Ay may be artificial because we have ignored 
CO photo-dissociation. The water ice abundance is not shown because it is almost constant 
2 X 10~^ throughout the plot. 
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Fig. 4. — Gas-phase H2O column density measured perpendicular to the disk midplane as 
function of the radius. In the very inner regions the H2O column density decreases because 
the gas densities at Aj, ~ 1 are higher, leading to a greater proportion of H2O that is 
frozen out onto grains as opposed to photo-desorbed. Once the grains are hotter than 100 K 
{R < 30AU), water will not freeze-out and gas phase abundances will rise again (not shown). 



